BEREZIN INTEGRALS AND POISSON PROCESSES 
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Abstract. We show that the calculation of Berezin integrals over anticommuting vari- 
ables can be reduced to the evaluation of expectations of functionals of Poisson processes 
via an appropriate Feynman-Kac formula. In this way the tools of ordinary analysis can be 
applied to Berezin integrals and, as an example, we prove a simple upper bound. Possible 
applications of our results are briefly mentioned. 
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1. Introduction During the last decades the functional integral has become the stan- 
dard approach to the quantization of systems with inflnitely many degrees of freedom. 
In typical cases like QED and QCD which involve both bosons and fermions the integral 
has to deal with anticommuting variables belonging to a Grassmann algebra and a lucid 
exposition of the rules of integration over these variables can be found in [Bl], [FS]. 

With the discovery of supersymmetry (SUSY) anticommutative integration has received 
further impetus and has been applied in different areas of physics and mathematics. SUSY 
has been flrst introduced in particle physics to express a possible fundamental symmetry 
between bosons and fermions and then has found several applications as a formal tool in 
the theory of complex systems like heavy nuclei and more recently disordered or chaotic 
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mesoscopic systems [VWZ, F, Z]. In mathematics it plays a role in various approaches to 
index theorems and other topics of differential geometry [BZ]. 

In physical applications of SUSY one introduces fields \1/ which have two types of com- 
ponents, boson-like <fi, ...,(f2n which are just ordinary numerical variables and fermi-like 

••••,^2n where the ^-s are anticommuting Grassmann variables. Supersymmetry trans- 
formations mix boson and fermion components. 

One is usually interested in the calculation of formal expressions like 

/ d[<i>mi[^ii[^jexp{-sm 

called correlation functions. <S'(^') is a functional invariant under supersymmetry transfor- 
mations. 

Integrals over anticommuting variables were introduced by Berezin also in statistical 
mechanics to represent the partition function of the planar Ising model [B2] and the gen- 
erating function of related combinatorial problems [B3]. See also [RZ]. 

The rules of calculation of these integrals, known as Berezin integrals, are recalled in 
the next section. In spite of the fact that this formalism allows compact and very powerful 
manipulations of the expressions of interest and in some cases exact calculations, it has 
the drawback that the usual tools of analysis are not applicable to the anticommuting 
variables. In particular one cannot easily find bounds for anticommutative integrals. 

The purpose of this paper is to show that the any Berezin integral can be represented 
in terms of the expectations of appropriate functionals of Poisson processes. On the basis 
of this representation ordinary analysis can be used and in particular upper bounds can 
be obtained. Furthermore correlation functions as above can be expressed entirely as 
expectations over ordinary stochastic variables. 

The starting point of our analysis is a generalized Feynman-Kac formula developed in 
the early eighties to express the solutions of the imaginary time Pauli equation [DJLS]. 
The equation is 

dti^t = -^(-^V - A{x)f,Pt - V{x)iPt + ^ll{x)-aiPt (1.1) 

A and H are the vector potential and the magnetic field respectively, a denotes the Pauli 
matrices in the usual representation. In [DJLS] we proved that the initial value problem 
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is solved by 

ipt{x, a) — 'E[ipo(x + wt, exp — / V{x + Wr)dT — i / A{x + Wr)-dwr+ 

Jo Jo 

(1.2) 

H,{x + Wr){-)''-adT + log[^{H,{x + Wr)-i{-f-aHy{x + Wr))]dNr)] 

a is a dichotomic variable which can take the values ±1. The expectation is taken with 
respect to the Wiener process Wt and to the Poisson process A^^. For the understanding of 
this formula we have to explain the meaning of the stochastic integral j dNf. A Poisson 
process is a jump process characterized by the following probabilities 

P(A^t+At -N, = k) = ^e-^^ 

Its trajectories are therefore piecewise constant increasing functions and we shall assume 
that at each jump they are continuous from the left. The stochastic integral is just an 
ordinary Stieltjes integral 



/ /(T)c/iV, = 5]/(r,) 

^0 I 



where are random jump times in the interval [0, t) which are distributed exponentially, 
that is P(r < t) = 1 - e"*. 

The interesting property of formula (1.2) is that by letting the spinor indices to become a 
stochastic process the Pauli matrices have disappeared from the expression of the evolution 
operator and their algebra is taken into account completely by the expectation over the 
jump process. The power of the approach was demonstrated by proving a non trivial 
paramagnetic inequality which shows that in three dimensions the evolution is bounded 
above by the evolution in a magnetic field which lies in a plane and whose components are 
simply related to the original magnetic field. This is easily seen by taking absolute values 
and implies for the ground states 

£;o(0, 0, 0; {Hf + Hlf/\ 0, H^) < Eo{A^, A^, A^, H^, H^, H^) 



For a recent application of the [D JLS] approach see [ER] . 
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Since the Pauli matrices are objects which belong to a Chfford algebra the above find- 
ings suggested a possible connection between calculus with Poisson processes and calculus 
with anticommuting variables. It is the main purpose of this paper to implement such 
a connection. To illustrate how the connection comes about in Sec. 3 we take up again 
the case of Pauli type equations and we observe that they can be interpreted as evolution 
equations over a Grassmann algebra. This type of evolution was considered for instance 
by Berezin and Marinov [BM] , [M] . Solutions of evolution equations over Grassmann alge- 
bras can be expressed as Berezin integrals which represent the convolution of the kernel of 
the evolution operator with the initial condition. A straightforward comparison with the 
solution given in [DJLS] provides the identification of the anticommutative integrals with 
the appropriate Poisson expectations. 

In order to develop the theory in a systematic way for an arbitrary but finite number of 
anticommutative variables in Sec. 4 we introduce a representation of Grassmann algebras 
in the space of functions of dichotomic variables that we call the cr-representation. This 
space was used by Wigner in his book on group theory [W] to find the representations of 
the symmetric group connected with the exclusion principle. 

In Sec. 5 we develop the necessary theory of semigroups associated with Poisson pro- 
cesses and in Sec. 6 we make the identification of Berezin integrals with appropriate 
expectations. We also derive a general inequality. 

In Sec. 7 we discuss in particular the Gaussian anticommutative integrals due to their 
importance in physical applications. 

We conclude this introduction with some comments on possible interesting applications 
of the results obtained in this paper. The semigroups associated to Poisson processes 
encompass all Hamiltonian semigroups {exp(— ti?)}^>Q for k interacting |-spins in an 
external magnetic field as, for instance, a Heisenberg ferromagnet or, by a slight change 
of language, any Hamiltonian semigroup for models describing interacting fermions on a 
finite lattice. Our representation can therefore be used both for theoretical or simulation 
studies of the statistical mechanics of such systems. 

A particularly interesting case to which our approach can be applied is the calculation 
of the Dirac propagator on a lattice, an important problem in the study of QCD on the 
lattice. Presumably in this way it is possible to obtain a simplification of the methods 
used at present in simulations. See the remark at the end of section 7. 

More generally in all cases where SUSY is relevant our approach may be a useful tool. 
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We believe that our results have also an independent interest in so far as they establish 
a direct connection between algebraic objects as those represented by Berezin integrals 
and analytic expressions. 

As a final remark we observe that the Wiener process and the Poisson process are both 
Levy processes and are actually the two limit cases of the Levy-Khinchin formula. We find 
therefore quite satisfying their correspondence at the euclidean level with the two basic 
types of particles in nature, bosons and fcrmions. Then the question naturally arises: 
do the other processes described by the Levy-Khinchin formula have any relevance for 
physics? 

2. Analysis on Grassmann Algebras 

We open this Section by a short review of standard definitions and results and refer 
reader to [B] , [D] , [FS] for more detailed information. 

By G{k) we denote the Grassmann algebra over C generated by its identity 1^ and a 
family {^i, ... of generators which obey the following commutation relations: 

CiCj = -^jCi; (2.1) 

For future reference, G^{k) ~ G{nk) will be the Grassmann algebra generated by 1^ 
and {^1, . . . , . . . , ^2 ^ . . . , . . . , }, and by convention, (k) = G{k). 

The collection {!,..., A; } of labels is any non - empty finite set endowed with a total 
ordering. 

Elements of G(k) of the form ^i^^^i^ • . . . • are called monomials; we will use the set of 
ordered multiindices Mj. = {fj, = {/ii, . . . , : 1 < //i < /X2 < • • • < A*n < k}, and write: 

for fj, = {/ii, ij,2, ■ ■ ■ ,fJ"n)- As a linear space G(k) has dimension 2^, and each element 
F(^) e G(k) can be represented in a unique way as a polynomial with complex coefficients: 

k 

m = m = fo-ik + Y. E /n,...,^?^ ■■■■■^^.= E ■ e^ (2.2) 

r=l l<ii<...<ij.<k iJ.€Mn 

where G C, and therefore G{k) is naturally graded. It is advantageous thinking of F(^) 
as a "function" of the Grassmann "variables" ^i, ... , that is to say of the "Fermi field" 



Analysis over G{k) is based upon left derivatives ^ and Berezin integral, 

which are defined as follows: 

k 

(—1) SfJ^ji ■ • . . . • . . . • , 

J=l 

where the sign / over the generator means that it is omitted. 

To define an integral we introduce symbols d^i, . . . ,d^k satisfying the following com- 
mutation relations: 

where {a,b} = a ■ b + b ■ a, and define "basic" integrals: 

B rB 

dii = 0, / i^-ii d^i = ; 

if fij 7^ i for all j in fi. We extend the integral on G(k) by linearity and call it the Berezin 
integral. In general 

FiOd^k • . . . • rfCi =' f = fi,...,k e C , (2.3) 

and observe that up to a multiplicative constant, is uniquely defined as the only linear form 
over G(k) alternating under permutations of the Grassmannian variables. It transforms 
[B] as 

FiR^Vk^ = (detR) ■ J Fi^Vk^, (2.4) 

under the linear substitution — > Yl'j=i ^ij^j- 

A specially important case is provided by "Gaussian integrals" : 

/B 1 ^ 

exp{- ^ Aij^i^jjVk^; (2.5) 
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where one can always assume that A = (Aij) is an antisymmetric matrix, otherwise it could 
be replaced by 2~^(A - A-^). "Gaussian integrals" vanish for odd k while, by exploiting 
(2.4), we get 

/ exp{ Yl AhkM'D^k^ = PfA-, (2.6) 

l<h<k<2k 

The Pfaffian Pf A of the triangular array {A-hk} i<h<k<2k defined by : 

PfA = J2{-lrA^,J, •...•A,,,,; (2.7) 

where the sum is taken over all {2k)\/2^k\ ways of pairing of the elements of the set 
{1, 2, . . . , 2k} where (— l)'^ is the parity of the permutation tt = (iiji, . . . , ikjk)- 

Let us consider the operator of left derivative ^ and the operator |j of left multiplica- 
tion by the element ^j, both acting on G(k) (see [B]). In the sequel we shall omit the hat 
if confusion does not arise. 

We recall that all linear operators acting on G{k) belong to the Clifi'ord (or Spinor) 
algebra C(2k) generated by its identity l2k and the operators . . . , and ^i, . . . , ^fe, 
which satisfy the following commutation relations: 

The algebra C(2k) is isomorphic to the CAR algebra [BR] over a 2^^ - dimensional Hilbert 
space since operators -A- and f j might be interpreted as annihilation and creation operators 
for a Fermi system with k degrees of freedom. 

We end the Section by defining kernels of operators acting on G(k). We recall that 
to each L G C(2k) corresponds a unique element Ker(L)(^,^') of the Grassmann algebra 
G(2k), generated by ^i, . . . , ^fe; and . . . , such that 

(LF)(0 = / Ker{m,e)-F{e)V,e, 
Jb 

3. Evolution on G{k) 
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Let us consider the evolution given by the equation: 

f = L/,; (3.1) 

where ft E G{k) and L e C(2k). 

In this section first we solve the equation (3.1) by constructing a kernel for the operator 
exp(tL) using standard tools of Grassmannian analysis and writing the solution of (3.1) as 
an element of G{k) with coefficients given by certain Berezin integrals, which we compute 
explicitely, and second, we show that these integrals could be represented as expectations 
with respect to a properly chosen family of standard Poisson processes. 

We introduce some necessary constructions related to the formal description of contin- 
uous time evolution on G{k). In order to do that it is convenient to embed G(k) into an 
extended Grassmann algebra Goo(^), and we will proceed in the following way (cf [ MIS]). 
Consider three Grassmann algebras: 

the Grassmann algebra Tr{k), t = (ti, . . . , tm), generated by: 

1, 6(ti), . . . ,^fc(ti); pi{ti), . . . ,pk{ti)\ i = 1, . . . ,m; 

the Grassmann algebra V^{k) generated by: 

1, 6('S), . . . ,^fc(s); Pi(s), . . . ,Pfe(s); s<t; 

and the Grassmann algebra V^oik) generated by: 

1, - Ak{s); pi(s),... ,Pk{s); s > 0; 

where 

Ci(s)-0(^) = -^jit) -^iis); Pi{s)-pj{t) = -pj{t)- pi{s), 

and 

Pi(«) -0(0 = -0(0 •Pi(«) "^hj, Vs,t>0. 

We consider exterior algebras: 

G,(A;) =Go(A;) A ^Ak); 

Gt]{k) = Go{k) /\Tt]{ky, 
Goo(A;) = Go{k) A r^oik); 

where Go (A;) = G(k). 

Remark, a) For the construction of Grassmann algebras with an infinite number of 
generators we refer also to [B], [R], [S]. b) The generators p's will be used only for "Fourier 
transform" - type expressions and play only an aiixilliary role here. 



9 



Example. The Pauli equation on G(l). 

Let us consider evolution on G(l), generated by 1 and ^i, given by (3.1) where 



L = hi ai + h2 (72 + hs 0-3 ; 



and we identify: 



so that the usual commutation rules of Pauli matrices are satisfied. 

As it will be proven in Theorem (3.1) the kernel of the operator exp{tL) can be written 
as the limit: 



Ker(e*^)(^,0= lim QT{^.0\ 



where 



Qr(e,o = f ■■■■■ f pn^,^t,Pt,e)T>mPtvut; 



m— times 



which is the kernel of the operator (1 + ;^L) and -P/^(^, ^t, pt, $,') G Go(l) A ^T{t,m) 
with T{t, m) = {t/m, 2t/m, . . . , [{m — l)t\/m, t), and has the following form: 



m 



m 



+ i/i2 



tl 
m 



tl 
m 



+ 



m 



m 



- P 



tl 
m 



m 



where 



T^mPt = dp 



m 



d^ 



dp{t); 



(m - l)t 



m 



m 
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and we put ^(0) = ^{t) = ^. In this way the solution of Equation (3.1) with initial data 
= /o-l + /i ■ /o, /i e C can be written: 

FtiO = (e'^^m) = f [Ker(e*^)(e,r)i^(e')] d^'; (3.2) 



and we have: 



where 



i^t(0 = /oWi + /iW-e, 



m = -ff [Ker(e*^)(e,^')i^(^') " ^] (3.3) 

= [Ker(e*^)(e,e')i^(e')] ^^e- (3.4) 

On the other hand from equation (1.2) specialized to the present simplified case we have 
that: 

fo{t) = e'E\f , 

2 

•exp(/ en{hi - ii-l)""- h2)dNr - [ hsi-l)''- dr))]; (3.5) 
Jo Jo 



and 



/i = e* E [/ i+(_i)ivt 

2 

•exp(/ en{hi+i{-lf-h2)dNr+ [ hsi-lf^dr)]; (3.6) 
where is a Poisson process with unit parameter. Comparing (3.3) and (3.5) we get: 



-J J [Ker(e*^)(^,Oi^(0^] d^ d^ = e'E[f^^^- 

• exp ( [ en{hi - h2)dNr - f hsi-l)^^ dr)] , (3.7) 

Jo Jo 

An analogous equation can be written for fi{t). 

In this way we obtain the equality between certain Berezin integrals and expectations 
with respect to the standard Poisson process. 
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To discuss the general case, let us consider the family of operators: 



which satisfy the commutation relations: 



j = 1,... ,n 



Further we will use the notations: letz^, x, /x G Mf. 



U = {Ui,... X ^ {Xi, . . . ,Xi}, /X = {/Xi, . . . 



such that 



we will write: 



i'r\ij, = i'r\x = sr\fi = 0; 



i=l J = l r=l 



Let us set: 



(3.8) 



{iy,x,ij.) 

where 7('^'»''M) e C{2k) and h^^^^^,,^ e C. 

Theorem 3.1. The kernel of the operator exp(tL) acting on G(k) is given as the limit: 



where 



Ker(e*^)(e,e')= lim QTi^.^) ^ Gi2k) ; 

m— »oo 

Qr(e,e') = 
• j](2e.p.+a)j]-(e^-zp^) 
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m 



my \ m 



Vp 



if k is even, and similar formula (see Appendix 1.) holds for k odd. 

Proof. The proof is rather simple and we postpone it to Appendix 1. 
In this way we can write solution of the equation (3.1) as the Berezin integral: 

with ff,(t) e C. 

Using Berezin integration rules we immediately get: 



= j [j [Ker(e*^)(e,e')i^o(e')]2^e'-e'^^}2^e; (3.9) 

where pf" e M is a complementary multiindex to /i, i.e. ji = {fxi,... , fin) and fx'^ = 
(//5,... , /i^,)> are such that {/ii,... flM^--- , ju^/} = and {/ii, . . . U M,--- 1 1^ 
{1,2,..., k}, and we assume here that $^ = {1, 2, k}. 

Let us fix a total ordering ^ on M (for example lexicographic order). It induces a total 
ordering on the set of monomials in the following way: 

if//^//'; 

let us rename all monomials with respect to the order ^ by . . . , which form the 

basis of G{k) as a 2'^-dimensional linear space. 

In this basis the matrix elements of the operator ^('^'^'/*) = 7j^. • 11^=1 (7a;j 7xj ) • 

Y[T=i explicitely computable, so equation (3.1) could be rewritten: 

2 ^ 

^ = E^«./5//3(^); (3-10) 

/3 = 1 
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where the coeflBcients -C/a,^ = L^^pihi, /i2, . . . , /i2fc) £ C. 

With an obvious definition of $(q;,/3) and '^{oi), equation (3.10) can be rewritten 



at 



2"- 1 



exp[$(a,/?)] • f^^f,{t) + - /alt); 



(3.11) 



with the initial condition 



/a(0) = /a; 

and where the sign © stands for the sum modufo 2^. 
The solution of the linear system (3.11) is given by 



Ut) 



/«® AT, • exp( *(a © iVt e iVr) + 



2'=-! 

+ ^ / $(a © iVi © iV^ © /3, 13) dN^ 



(3.12) 



where AT^ = E;s="i P^t is the sum of 2'^ — 1 independent Poisson processes and a © 
Nf, a ® Nt Q Nt-, and a ® Q N^- Q P stand for sums and differences modulo 2^^. (For 
complete proof see [DJLS].) 

Now comparing (3.9) and (3.12) 



fa®N,- ^^p(^J^ ® Nt e Nr) dr + 



,(2'=-l)-t 



+ J2 / Ha ® Nt e Nr e P, P)dN^ 
= 1 



(3.13) 



Sums and differences modulo 2^ are slightly unconfortable to handle. In the next three 
sections we shall reformulate the theory in a space of functions of dichotomic variables. 
This allows the construction of a systematic formalism suitable for both theoretical and 
numerical analysis. 
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4. Representation of Grassmann Algebras in the Space of Functions of Di- 
chotomic Variables (a - Representation) 

In this Section we discuss a linear bijection between Grassmann algebras and spaces of 
functions of dichotomic variables. 

Let Z2 be {—1, 1} with its natural Abelian group structure and Z2 be the direct 
product of k copies of Z2 which is finite commutative group with the unit element = 

Let us define Tik as the linear space of all "wave functions" x(-) : 22^^ — C of /c di- 
chotomic variables a"i, . . . , Ufe. It becomes - dimensional Hilbert space when equipped 
with the inner product: 

<Xi,X2>k= Yl Xi(f^)x2(f^); (4-1) 

and could be interpreted as the space of (pure) states for a Heisenberg ferromagnet in a 
finite box. 

All monomials of the Grassmann algebra G{k) can be indexed by elements of Z2 in 
the following way: 

a ^ = ^'^ = Ci^ •••4^, ^ = (f^i, ••• ,^fe) eZ2^^ (4.2) 
and let us define now a map X: 7ik —>■ G{k) by the formula: 

^(X) ^ i^x(e) = E x(^)-e'^ e Gik); (4.3) 

where the "wave function" x(.) G Hk- 

It is easy to see that the map X is injective, i.e. I{Xi ) 7^ ^(X2 ) if Xi 7^ X2 • Moreover, 
the map J : G{k) ^ Hk defined by: 

J(F(0) = X.(^) = Afc(a) f^'^-FiOm, (4.4) 
where Afc(.) : Zj ^ Z2 is given by 
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is the inverse of J : X = J~^, and, clearly, it is surjective. By this we obtain a linear bijec- 
tion between Grassmann algebra G(fc) and the space of functions of dichotomic variables 
Tifc, which we call a - representation. 

Next we turn to C{2k). To each A e C(2/c) corresponds a linear operator A : Hk 
7ik given by formula: 

AFiO = Yl (^X)H-^'^. (4.6) 



For any A e C{2k), which is a linear combination of normally ordered products ^ 2' 



l-?7 

(j^) ^ , where e, ry € Z2 , we find its image A by computing the image of operators 
^~2- ■ (A)-2^ which we denote by a*~2^ a~2^ . 
Proposition 4.1 For all e, ri e we have: 



where 



Ck{e,a) = n( 
1=1 ^ 



' 1 + ei I — ei 1 — (Ti 

— 7, \ Ti ■ Ti •€!••• ei-i • (Ti • • • (J/_i I ; 



M'n,<^) = n ( 

1=1 ^ 



7; 1 Ti Vl--- Vl-1 • C^l • • • CTl-l 



Nk{e,r),a) = Ak{r),a) • Cfc(e, r/cr) ; 

and where ecr = eiUi, . . . , ekCTk- 

Proof. Immediately follows from the fact that 

1-e , 1-17 , 1-eiT I S \ ^ . l-g , 1-T)g- 



= C'(e,e(7)-e^; — J = ^(r/, 77(7) • ; 

which could be checked by induction on k. 

l-e; 

Corollary. IfF(^) = Yliee'L^'' fi^)^~^ ^ G(/c), the image F (a*) of the linear operator 
F(0 e C(2/c) is 

(F(a*)x)H= E /(e)Cfc(e,a)x(ea), (4.8) 
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and, in particular, when F{^) = Ei<i<j<2fc ^ij^i^j 

^<i<j<2k 

= ^ ^ij — — n •••,-f^j,---,f^2fe), (4.9) 

{i,j)eT{A) i=i+i 

where r(^) = {(z, j), 1 < z < j < 2A; : }. 

Let ^ be a hnear operator acting on Ti^ and A{., .) its matrix, given by 

{Ax){a)= Yl M<r,<T')x{<T'). (4.10) 

For instance, from Proposition 2.1 we get that the matrix element of the operator a*~^a^^ 
is given by 

A{(T,(t') = Nk{e,r],(7)6ea,r,a' ■ 

Now from the formula (4.4) we obtain: 

p ^'^Ai^Vf.i = Ak{a)A{a,a'), e C{2k) (E.l) 

which relates the Berezin integral A^^'^V^^ to the matrix A{a, a'). This is our 

first basic formula. 

We end the Section by some remarks about kernels of ^ e C{2k). We recall that to 
each A e C{2k) corresponds a unique Ker(74)(^, ^') e G^(/c), the kernel of the linear 
operator A, such that 

AF{0 = f Ker{A){^,OF{i')T^k^' . VF(0 G G(/c). (4.12) 

By exploiting Lemma 2.1, it is easy to see that for all positive integer k the kernel 
Ker(l2fe)(C) CO of the unit element I2K of C{2k) is 



Ker(l2fc)(C,0 = E A,(a)C^C'^, 



(4.13) 
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and more generally: 

Ker(A)(C,0= E Afc(a')A(a,aOe^e"^, (4-14) 

<J,<7 fc.ti2 

which, together with 

A{a,a') ^ Ak{a) jj ^'^KeT{A){i,i') i di\- ■ ■ di\dik ■ ■ ■ di^ , (4.15) 

relates the kernel Ker(A)(^, of A and the matrix A(cr, cr'). 
5. Semigroups Associated to Poisson Processes 

We now turn to semigroups of linear operators acting on the hilbert space T-Lk- We 
shall give a probabilistic representation of the semigroup {exptLj^^Q generated by a non - 
trivial linear operator L : Hk Ti-k- It specializes the more general formulas introduced 
in [DJLS] to which we address the reader, nevertheless the Section will be self - contained. 

We start with a special representation of L G L{TCk, 'Hk)- For all e = (ei, . . . , 6^) G 
2i2^, and different from the identity : e = (1, . . . , 1), let D'^ be the self - adjoint 

difference operator: 

(D^xK^) = - x(^); (5.1) 

which annihilates constants. 

We observe that each L e L{Ti,k-, T^k) admits the representation 

L = 5^A,(.)D^ - F(.)l, 

where the functions Ae(.), V{.) : 7i2^ — > C are related to the matrix !/(.,.) of the 
operator L by formulas: 

i) K{a) = L{a,e<7); ii) V{a) = - J] L{a,a'). 

Let r(L) be the collection of e 7^ such that Xe{-) is not identically vanishing and 
|r(L)| be its cardinality. If r(L) ^ we call the operator L a difference operator. 
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Now let {N[}^_^^^ be a given collection of {2^ — 1) independent Poisson processes of 
unit parameter which we assume to be left - continuous. 

Theorem 5.1 (Probabilistic representation of semigroups) 
Let L = Sg£r(L) '^e(-)-^'^ ~ ^(O^ be a difference operator. Then 

= e*I^WlEfx(cx(-l)^*)exp{ ( [ A,(a(-l)^OdiV| - 

where iV, = (Ar^, . . . , N^) with iV] = Eeer(L) ^^I- 

Remark By convention exp J^^ /n Ae(cr(— l)'^^)(iAr| vanishes if Ae(cr(— 1)''^^) = for 
some < s < t such that Nl_^_ — iV| 7^ 0. We observe that it doesn't depend upon the 
choice of the branch of the logarithm. 

Proof. We follow the strategy explained in [DJLS]. 

i) The r.h.s. defines a semigroup {P*}^>q of linear operators on Tix by the Markov property 
of Poisson process. In order to complete the proof we must show that the infinitesimal 
generator of {P*}^>q coincides with L. 

ii) Since each wave function x(.) is a linear superposition of characters Xn(c) = a"^ = 
cr"^ . . •cr^'", n = {ni, . . . ,nk) G (Z^^)^^ , Z2 = {0, 1}, we can only consider the case of 
x(-) = Xn(-) for some n e (73)^'^. Let ^^'^ be the random variable 

eer(L)'^[0'*) -^0 

where 

be,n{^) = i^^—^m +lnXe{cr); 
1=1 

mia) = - E K{a) - V{a) . 
eer(L) 

By definition, 

{P\n)ia) =Xn(^)e*l^(^)lE(exper). 
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iii) The stochastic differential dexp^^'^ of the process t e [0, +00) — > exp^^'" can be 
exphcitly evaluated as explained in [DJLS]. It turns out that 

c^exp^^''''' = 

= (exper")[m(a(-l)^*)dt + ^ {e^^--^''^-^^'''^ - 1) dN^] = 

eGr(L) 



(exp^r)[-^(^(-l)'^*)^^ + E A,(a(-l)^*)(e-cZAr,^ - dt) - dN, 

eeT{L) eer(L) 



By taking the expectation of dexp^^'" for t = 0, since E{dNl) = dt and {D^Xn){'^) = 
(€"■ — l)xn){(^), it follows that ■^(P*Xn)((^)\t=o = iLxn)i(^)- Therefore the infinitesimal 
generator of the semigroup {-P*}t>o exactly the difference operator L. 

Example As an elementary illustration of Theorem 5.1, let k = 1 and {Dx){cr) = 
x(-f^) - X(f^), then 

(e*^x)(^) =E(x(a(-l)^*)), (5.2) 
where Nf is the Poisson process with unit parameter. Indeed 

00 —tj.n 

ft* 

n=0 

= e-Hx(f^)cosht + x(-(^)sinht} = (e*^x)(f^)- 
From Theorem 5.1 we obtain the matrix elements of the operator e*^: 

- j\,{a{-lf^)ds^ - j\{a{-lf^)ds}^ , (£;.2) 

which is the second important equality and will be used in the next Section evaluating 
Berezin integrals. 
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6. Berezin Integrals and Poisson Processes in the a - Representation 

We shall consider Berezin integrals of the form: 



J rexp(-5(0m = y" er^2^---Cexp(-5(0m, (6.1) 

where n = (ni, . . . , Uk) G (Za)^'^ and S{$,) 7^ c • 1^ is a non - trivial element of G{k). It 
could be interpreted as the "action" for the (Euclidean) "Fermi field" {^1, . . . , ^k} in which 
case (6.1) would provide all (unnormalized) "correlation functions" or (Euclidean) "Green 
functions" of the field. In particular, for n = (0, 0, ...,0), (6.1) gives the "partition 
function" : 

Z[S] = j"" eM-SmT^k^- (6.2) 

Remark. We could easily consider more general integrals by taking for S an element of 
C(2k). Integrals of the form (6.1) cover a large class of physical applications. 
Let s(.) : Z2 ^ C be defined as 

s{e) = Afc(e) r^'^SiOV.C. (6.3) 



with Ak{e) as in (4.5). 

Theorem 6.1 (Berezin integrals as Poisson averages) For each /c = 1, 2, . . . , each non 
trivial S{^) G Gfc and for all n = (ni, . . ., Uk) e {7i\)^^ 

-B 



exp{ J] / ln{-s{e)Cu{e,-{-lf^^^))dNl}), 



X 

eer(5)'^[°'^) 
where s(e) is given by (6.3). 
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In particular 



Z[S] = e^\ns)\-s{e,))^fY[ 

^1=1 



k_ ^ _ Ni 
X 



Proof. Using formula (4.4) we get: 



exp{ J2 I ln{-s{e)C^{e-{-lf^))dNl}] 
eer(S)>^[0'i) ^ 



xexp{ V / ln{-s{e)Cu{e,-{-lf^))dNl}]. (6.5) 



and, by hypothesis, the subset T{S) = {e e 712^,6 ^ ek : s(e) ^ 0} is non empty. To the 
operator — <S'(^) e C{2k) corresponds the image L = —S{a*) e L{Hk,'Hk) which is the 
difference operator: 

L = E - V{.)1, (6.7) 

eGr(5) 

with 

X,{a) = -s{e)Ck{e,a), (6.8) 
V{a) = E <^)Ck{e,^) = s{ek) - E Xeia). (6.9) 

From equations (E.l), (E.2) and equality 

-V{a) - E A,(a) = -s(efc), 
eer(S) 



follows that 



Jo ^-n/c'\JO 



€Gr(S) 



Cfc(e, 0") is defined in Proposition 4.1 and Ng, in Theorem 5.1. 
We now make an important remark. 
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Remark . In equation (E.3) only the trajectories of the Poisson process with zero or one 
jump contribute to the expectation. In fact as soon as one of the factors (1 + (—1)^* 
in Ck vanishes the stochastic integral equals — cxo. On the other hand A''^'^ is a sum of 
independent processes and the event of two processes jumping at the same instant has 
zero probability. 

We now take advantage of the fact that the calculation of a Berezin integral has been 
reduced to an ordinary integral and we derive simple estimates. From Theorem (6.1) we 
get 



< 



< g(|r(5)|-fles(efc))g/ -Q 



1 _ 



X 



X 



exp{ V / /n(|s(e)| 
eer(5)-^[0'i) 



Cfe(e,-(-l)^^+^) 



(6.10) 



k 

x{Ni=o,i} n 

1=1 



N[+n, 



n W« 

eeris) 



,i-nesie.))l_ Y: (-1) 
pi,....pfc=0,l 



eer(5) 



since |Cfe(e, cr)| < 1 and J^Qi^dNg = A^;^; x is the characteristic function of the event 
indicated. 

7. Gaussian Berezin Integrals 

Gaussian Berezin integrals are a particular case when the "action" S is bilinear in the 
"Fermi field". Let G'^{k) be the Grassmann algebra generated by {^i, ■ ■ ■ ■,f]k'i'niT • ■ iVk} 
and {B)ij any k x k matrix. Let us consider the integral 



/B \ ^ 

Vi' ■■■ rjk'rfi ■ ■ ■ ^'"expl-- BijfjirjjjVkVT^kV , (7.1) 
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where z^/jZ^; e Z2 = {0; !}• Making the substitution ^1 = f]i,...,^k — Vk', Ck+i = 
?7i,---,6fc we get 

k 2k 

where (A^s) is the 2k x 2 A; antisymmetric matrix defined by 

2~^S^(s_/j) when r = 1, . . . , A; and s = k + 1, . . . , 2k; 
^rs = i —'^~^Bg(^r-k) when r = A; + 1, . . . , 2A; and s = 1, . . . , A;; (7.2) 
otherwise. 

Therefore the Gaussian integral (7.1) can be written in the foUowing form: 
,B 1 k 

Y. 



/B k 



/B 2fc 
^r---aexp(-- Yl Ars^rL)'^2k^ ^ (7.3) 

and from now on we use the last form. Let r{A) be the set: 

r(^) = {(r, s), l<r<s<2k : 7^ 0} , (7.4) 

which we suppose not to be empty. From (4.9) we get that the generator L of the corre- 
sponding semigroup is given by 

(LV^)(a) = 

1 — (T 1 — a 

= Y — — ^ JJcriV'(a-i,...,-c7r, ...,-crs,...,a-2fe), (7.5) 

{r,s)er{A) l=r 

and therefore 



2k 



r,s=l 
2A: . / .NAr{+ni 



A.(-(-i)>i-(^)iE(n^^^^ — n (^^^)"^ X 

^1=1 {r,s)eT{A) 



{EA) 
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X exp{ ^ / In 
(r,s)er(A)-^[0'i) 



l=r 



where {N^ ' } s)er{A) ^ family of independent Poisson processes with unit parameter 
and 



T{m,l) 



for all 1< / < 2/c. 



(z,m)er(A) (m,Oer(A) 
Example. Let k = 1, n= (0,0) and (^rs) = i(^2 (Pauli matrix), then 

/B pB 



eE 



l-(-l) 



exp / In 

'[0,1) 



-(-1) 



(7.6) 



This equality can be checked directly: indeed the random variable under expectation 
vanishes unless the Poisson process has exactly one jump for some < s < 1, in which 
case A^i = 1, and which happens with probability . 

Remark . Suppose we want to use our representation of Gaussian Berezin integrals 
in a numerical simulation. Let 2k be the total number of r] and fj. The generic case 
will require the generation of k{2k — 1) independent Poisson processes. For a Dirac field 
on a finite d-dimensional lattice A the maximum number of Poisson processes involved is 
L|A|(2L|A| — 1)) where L is the number of components of the field. This is a large number. 
However if we want to calculate for example the propagator of the free field we need much 
less. In fact the matrix Bij which is the lattice version of the differential operator ^-\-J^ — M 
couples only nearest neighbours. Therefore the required number of independent Poisson 
processes is reduced to roughly dL\h\ since each site has 2d nearest neighbours. 
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Proposition A.l. The kernel of the operator ^('^^'^'i^^ given by (3.8) and acting on 
G{k) is given by the foUowing formulas: 

/B n £ 



m ^ k 

if k is even, and 

/B n 



• n i^^M. - Pn) ■ { - 5^Pr(er + C)}dp; iA.2) 

j=l r=l 

if k is odd. 

The proof is straightforward, and we omit it. 

Remark. Here we expressed elements of the algebra Go{k) in so called "Fourier- 
transform" form. In fact we just multiply initial expression by p-monomials and then 
integrate over p. 

Let us set: 

(i/,a;,/[i) 

where 'j^'''^'^^^ g C(2A;), and are of the form (3.8), and hf^y^^.^^^ G C. 
After computations using (A.l) and (A.2) we get: 

Ker(e^ ^)(e, O = Ker(l - O + o(-) = 

n n 

/B 2 ^ _— _ _— _ 

®^Pi~^ XI h^,x,^l)W{^l^ + ^Pv)W{'i^xPx + ^'^)■ 

{u,X,fJ,) V X 

^ n ^ 

• n ~ -^pAir + 0}dpn • . . . • rfpi + o(-), 

[X r=l 
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for k even and analogous formula holds for k being odd. 

Remark. Since n is fixed we always understand convergence as a pointwise convergence 
in a finite dimensional linear space. 

Finally applying Trotter's formula we have: 



/B pB 
... Ker(e^^) 
J 

•Ker(e'" ) 



{m-l)t 



m 



(m-l)t \ ^({m-2)t 



m 



m 



. . .Ker(e"^-^) 



m 



{m -l)t 



m 



"B /-B f / j ■ t 



•j](2e,p,+ii)-n-(^M-^PM) 



i=i ''=1 



m 



f-) • • • • • dp{t)di'di 

= hm r Qn^,e)f{e)d^'; 

m— >oo J 

and from which we get Theorem 3.1. 



dm 



{A.3) 
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